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Axilp sequence folowng Ser 306 and occurs wrtftin 
the domain of A)dlp that shows homology with hIDE 
(74). To delete the complete STE23 seguence and 
create me ste23X:URA3 mutation, polymerase chain 
reaction (PCR) pnmers (5 ' -TCGGAAGACCTCAT- 
TCTTGCTCATTTTGATATTGCTC* TGTAGATTG- 
TACTGAGAGTGCAC-3*; and S'-GCTACAAACAGC- 
GTCGACTTGAATGCCCCGACATCTTCGACTGT. 
GCGGTATTTCACACCG-3') were used to ampfify 
the URA3 seguence of pRS3i6. and the reaction 
product was tran sformed into yeast tor one-step gene 
replacement [R Rotnstem, Methods Enzymol. 194, 
281 (1 991 )]. To create the axil A::LEU2 mutation con- 
tained on pi 14, a 5.0-kb Sal I fragment from oA&7 
was ctoned into PUC19, and an intemal 4.0-kt> Hpa 
^Xho I fragment was repiaced with a LBJ2 fragment 
To construct the ste23A::LEl!2 allele (a deletion cor- 
responding to 931 amino adds) earned on p153. a 
LBJ2 fragment was used to replace the 2.6-to Pml 
^-€d136 II fragment of STE23. which occurs within a 
6.2-kb Hind W-Bgi II genomic fragment carried on 
pSP72 (Promega). To create YEpMFAl, a 1.6-kto 
Bam HI fragment containing MFA1. from pKK16 [K. 
Kuchler, R. E. Sterne, J. Thomer, BMBO J. 8, 3973 - 
(1 989)]. was ligated into the Bam HI site of YEd351 [J. 
E. H3l, A. M. Myers. T. J. Koemer, A. Tzagotoff, Yeast 
2.163(19865]. 

24. J. Chant and I. Herskowte, Cell 65, 1203(1991). 

25. B. W. Matthews. Acc. Chem. Res. 21. 333 (1988). 

26. K. Kuchler, H. G. Dontman. J. Thomer. J. Cell So/. 
120. 1203 (1993); R. Kotting and C. P. Hoflenberg. 
BMBO J. 13, 3261 (1994); C. Berkower, 0. Loayza, 
S. Michaelis, Mol Biol. Cell 5, 1185(1994). 

27. A. Bender and J. R. Pringte. Proc. Natl. Acad. Set. 
USA 86. 9976 (1989); J. Chant. K. Corrado, J. R. 
Pringte. I. Herskowitz, Cell 65, 1213 (1991); S. 
Powers, E. Gonzales, T. Christensen. J. Cubert, D. 
Broek. ibid., p. 1225; H. O. Park. J. Chant, I. Her- 
skowitz. Nature 365, 269 (1993): J. Chant. Trends 

Genet 10, 328 (1994); and J. R. Pringle. J. 

Cell Biol 129. 751 (1995); J. Chant, M. Mischke, E. 
Mitchell. I. Herskowitz. J, R. Pringle, ibid., p. 767. 

28. G. F. Sprague Jr., Methods. Enzymol. 194. 77 
(1991). 

29. Single-letter abbreviations for the amino acid resi- 
dues are as follows: A. Ala; C, Cys; D. Asp; E. Giu; F, 
Phe; G. Gry; H. His; I. lie; K, Lys; L Leu; M. Met; N, 
Asn; P. Pro: Q. Gin: R, Arg: S, Ser; T. Thr; V. Val; W. 
Trp; and Y, Tyr. 

30. A W303 1A derivative. SY2625 (M47a ura3-T teu2*3. 
I12trpi-1 ade2- 1 canui00sst1& mfa2^:FUSl- tacZ 
rus3£L-:FVS I -H/S3), was the parent strain fa the mutant 
search. SY2625 denvatives for the mating assays, se- 
r-reted pheromone assays, and the pulse-chase exper- 
iments included the following strains: Y49 {ste22-1) t 
Y115 (mfe;A.;iH/2), Y142 iaxJ1 ::URA3), Y173 
iaxl11::LBJ2), Y220 iaxl1::URA3 ste23£L':URA3l. Y221 
iste23£L-:URA3), Y231 iaxl1L':L£U2 ste23A::LEU2}, 
and Y233 iste23t::LBJ2\. MATa denvatives of 
SY2625 included the following strains: Y199 
(SY2625 made MA 7a), Y278 \ste22-1). Y195 
(mfaW:LEU2), Y196 {axil £L"LBJ2), and Y197 
(axil :;URA3). The EG123 (/WW 7a Ieu2 ura3 trp1 can) 
hts4) genetic background was used to create a set of 
strains for analysis of bud site selection. EG 123 de- 
rivatives included the following strains: Y175 
(axt1k:LEU2). Y223 {axl1::URA3). Y234 (sfe23A.v 
LEU2), and Y272 {axJ1A::LEU2 ste23A::LEU2). 
MATa denvatives of EG123 included the following 
strains: Y214 (EG123 made MATa) and Y293 
(ax/7A:.l£U2). All strains were generated by means 
of standard genetic or molecular methods mvoMng 
the appropriate constructs (23). In particular, the ax/7 
ste23 double mutant strains were created by cross- 
ing of the appropriate MATa ste23 and MATa axil 
mutants, followed by sporutation of the resultant dip- 
loid and isolation of the double mutant from nonpa- 
rental di-type tetrads. Gene disruptions were con- 
firmed with either PCR or Southern (DNA) analysis. 

31. P129 is a YEp352 p. E HO, A. M. Myers, T. J. Ko- 
emer. A. Tzagotoff . Yeast 2, 1 63 (1 986)] piasmid con- 
taining a 5.5-kb Sal I fragment of DAX17. p15l was 
derived from pi 29 by insertion of a linker at me Bg) If 
site wrthn AXL 7 , which led to an rvframe insertion of 
the hemagglutinin (HA) epitope (DOYPYDVPOYA) (29) 
between amno acids 854 and 855 of the AXL / prod- 



uct. pC225 is a KS+ (Stratagene) piasmid comartng 
aO^^BamHl-SstlfragmentfrcfnpAXL7.SuPsb- 
tuoon mutations of the proposed active site of AxJtp 
were created wrm the use of pC225 and site-specrfic 
mutagenesis invoiving appropriate synthetic ofagonu- 
cteotides &dl-H6&A, 5 ' -GTGCTCACAAAGCGCT- 
GCCAAACCGGC-3'; axtUE7lA. 5'-AAGAATCAT- 
GTGCGCACAAAGGTGCGC-3 ' ; and axJ1~E7lD. 5'- 
AAGAATCATGTGATCACAAAGGTGCGC-31 The 
mutations were confirmed by seouence anafysis. Af- 
ter mutagenesis, the 0.4-kto Bam Hl-Msc I fragment 
from the mutagemzed pC225 piasmids was trans- 
ferred into pAXL 1 to create a set of pRS31 6 ptasmtis 
carrying different AXL 7 alleles. P124 Iaxf1-H68A), 
P130 {axfl-E7iA), and p132 <ax/I-E77D). Similarly, a 
set of HA-tagged atetes carried on YEp352 were cre- 
ated after replacement of the pl5i Bam HMvtsc I 
, to generate pl6i tax? 7 -£77 A), pi 62 (ax/7- 
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Quantitative Monitoring of Gene Expression 
Patterns with aJComplementary DNA Microarray 

Mark Schena,* Dari Shalon,*t Ronald W. Davis, 
Patrick O. Brown* 

A high-capacity system was developed to monitor the expression of many genes in 
parallel. Microarrays prepared by high-speed robotic printing of complementary DNAs on 
glass were used for quantitative expression measurements of the corresponding genes. 
Because of the small format and high density of the arrays, hybridization volumes of 2 
microliters could be used that enabled detection of rare transcripts in probe mixtures 
derived from 2 micrograms of total cellular messenger RNA. Differential expression 
measurements of 45 Arabidopsis genes were made by means of simultaneous, two-color 
fluorescence hybridization. 



The temporal, developmental, topographi- 
cal, histological, and physiological patterns 
in which a gene is expressed provide clues to 
its biological role. The large and expanding 
database of complementary DNA (cDNA) 
sequences from many organisms ( I ) presents 
the opportunity of defining these patterns at 
the level of the whole genome. 

For these studies, we used the small flow- 
ering plant Arabidopsis thaliana as a model 
organism. ATobidopsis possesses many ad- 
vantages for gene expression analysis, in- 
cluding the fact that it has the smallest 
genome of any higher eukaryote examined 
to date (2). Forty-five cloned Arabidopsis 
cDNAs (Table 1), including 14 complete 
sequences and 31 expressed sequence tags 
(ESTs), were used as gene-specific targets. 
We obtained the ESTs by selecting cDNA 
clones at random from an Arabidopsis 
cDNA library. Sequence analysis revealed 
that 28 of the 31 ESTs matched sequences 
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in the database (Table 1). Three additional 
cDNAs from other organisms served as con- 
trols in the experiments. 

The 48 cDNAs, averaging ~L0 kb t 
were amplified with the polymerase chain 
reaction (PCR) and deposited into indi- 
vidual wells of a 96-well microliter plate. 
Each sample was duplicated in two adja- 
cent wells to allow the reproducibility of 
the arraying and hybridization process to 
be tested. Samples from the microttter 
plate were printed onto glass microscope 
slides in an area measuring 3.5 mm by 5.5 
mm with the use of a high-speed arraying 
machine (3). The arrays were processed by 
chemical and heat treatment to attach the 
DNA sequences to the glass surface and 
denature them (3). Three arrays, printed 
in a single lot, were used for the experi- 
ments here. A single microtiter plate of 
PCR products provides sufficient material 
to print at least 500 arrays. 

Fluorescent probes were prepared from 
total Arabidopsis mRNA (4) by a single 
round of reverse transcription (5). The Ara- 
bidopsis mRNA was supplemented with hu- 
man acetylcholine receptor (AChR) mRNA 
at a dilution of 1 : 10,000 (w/w) before cDNA 
synthesis, to provide an intemal standard for 
calibration (5). The resulting fluorescently 
labeled cDNA mixture was hybridized to an 
array at high stringency (6) and scanned 
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with a laser (3). A high-sensitivity scan gave 
signals that saturated the detector at nearly 
all of the Arabidopsis target sites (Fig. 1A). 
Calibration relative to the AChR mRNA 
standard (Fig. 1A) established a sensitivity 
limit of - 1 : 50,000. No detectable hybridiza- 
tion was observed to either the rat glucocor- 
ticoid receptor (Fig. 1A) or the yeast TRP4 
(Fig. 1A) targets even at the highest scan- 
ning sensitivity. A moderate-sensitivity scan 



High sensitivity 
4 5 6 7 B 9 10 11 12 

•r. p o ^ ^ -3 *v 



of the same array allowed linear detection of 
the more abundant transcripts (Fig. IB). 
Quantitation of both scans revealed a range 
of expression levels spanning three orders of 
magnitude for the 45 genes tested (Table 2). 
RNA blots (7) for several genes (Fig. 2) 
corroborated the expression levels measured 
with the microarray to within a factor of 5 
(Table 2). 

Differential gene expression was investi- 
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gated with a simultaneous, two-color hv- 
bridttation scheme, which served to mini- 
raize experimental variation inherent in the 
comparison of independent hybridizations. 
Fluorescent probes were prepared from two 
mRNA sources with the use of reverse tran- 
scriptase in the presence of fluorescein- and 
lissamine-labeled nucleotide analogs, re- 
spectively (5). The two probes were then 
mixed together in equal proportions, hy- 
bridized to a single array, and scanned sep- 
arately for fluorescein and lissamine emis- 
sion after independent excitation of the two 
fluorophores (3). 

To test whether ovcrexpression of a sin- 
gle gene could be detected in a pool of total 
Arabidopsis mRNA, we used a microarray to 
analyze a transgenic line overexpressing the 
single transcription factor HAT4 (8). Fluo- 
rescent probes representing mRNA from 
wild-type and HAT4-transgenic plants were 
labeled with fluorescein and lissamine, re- 
spectively; the two probes were then mixed 
and hybridized to a single array. An intense 
hybridization signal was observed at the 
position of the HAT4 cDNA in the lissa- 
mine-specific scan (Fig. ID), but not in the 
fluorescein-specific scan of the same array 
(Fig. 1C). Calibration with AChR mRNA 
added to the fluorescein and lissamine 
cDNA synthesis reactions at dilutions of 
1:10,000 (Fig. 1C) and 1:100 (Fig. ID), 
respectively, revealed a 50-fold elevation of 
HAT4 mRNA in the transgenic line rela- 
tive to its abundance in wild-cype plants 
(Table 2). This magnitude of HAT4 over- 
expression matched that inferred from the 
Northern (RNA) analysis within a factor of 
2 (Fig. 2 and Table 2). Expression of all the 
other genes monitored on the array differed 
by less than a factor of 5 between HAT4- 
transgenic and wild-type plants (Fig 1, C 



Wild type HAT4 transgenic 



CAB! 



HAT4 



ROC1 



Rg 1 Gene expression monitored with the use of cDrM microarrays. Ruorescent scans represented in 
pseudocolor correspond to hybridization intensities. Color bars were calibrated from the signal obtained 
with the use of known concentrations of human AChR mRNA in independent experiments. Numbers and 
letters on the axes mark the position of each cDNA. (A) High-sensitivity fluorescein scan after hybndization 
with fluorescein-labeled cDNA derived from wild-type plants. (B) Same array as in (A) but scanned at 
moderate sensitivity. (C and D) A single array was probed with a 1 : 1 mixture of fluorescein-labeled cDNA 
from wild-type plants and lissamine-labeled cDNA from HAT4-transgenic plants. "The single array was 
then scanned successively to detect the fluorescein fluorescence corresponding to mRNA from wild-type 
plants (C) and the lissamine fluorescence corresponding to mRNA from HAT4-transgenic plants (D). (E 
and F) A single array was probed with a 1:1 mixture of fluorescein-labeted cDNA from root tissue and 
lissamine-labeied cDNA from leaf tissue. The single array was then scanned successively to detect the 
fluorescein fluorescence corresponding to mRNAs expressed in roots (E) and the lissamine fluorescence 
corresponding to mRNAs expressed in leaves (F). 
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AChR ^^^^^ 

20 2.0 0.2 
mRNA (ng) 

Rg. 2. Gene expression monitored with RNA 
(Northern) blot analysis. Designated amounts of 
mRNA from wild-type and HAT4-transgenic 
plants were spotted onto nylon membranes and 
probed with the cONAs indicated. Purified human 
AChR mRNA was used for calibration. 



and D, and Table 2). Hybridization of flu- 
orescein-labeled glucocorticoid receptor 
cDNA (Fig. 1C) and lissamine-labeled 
TRP4 cDNA (Fig. ID) verified the pres- 
ence of the negative control targets and the 
lack of optical cross talk between the two 
fluorophores. 

To explore a more complex alteration in 
expression patterns, we performed a second 
two-color hybridization experiment with 
fluorescein- and lissamine-labeled probes 
prepared from root and leaf mRNA, respec- 
tively. The scanning sensitivities for the 
two fluorophores were normalized by 
matching the signals resulting from AChR 



mRNA, which was added to both cDNA 
synthesis reactions at a dilution of 1:1000 
(Fig. 1 , E and F). A comparison of the scans 
revealed widespread differences in gene ex- 
pression between root and leaf tissue (Fig. 1, 
E and F). The mRNA from the light-regu- 
lated CABi gene was -500-fold more abun- 
dant in leaf (Fig. IF) than in toot tissue 
(Fig. IE). The expression of 26 other genes 
differed between root and leaf tissue by 
more than a factor of 5 (Fig. 1, E and F). 

The HAT4 -transgenic line we examined 
has elongated hypocotyls, early flowering, 
poor germination, and altered pigmentation 
(8). Although changes in expression were 



Table 1. Sequences contained on the cONA rnicroarray. Shown is the position, the known or putative 
function, and the accession number of each cDNA in the rnicroarray (fig. 1). All but three of the ESTs used 
in this study matched a sequence in the database. NADH, reduced form of nicotinamide adenine 
dinucieotide; ATPase, adenosine triphosphatase; GTP, guanosine triphosphate. 



Position 


cDNA 


a1.2 


AChR 


a3. 4 


EST3 


a5,6 


EST6 


a7, 8 


AACi 


a9, 10 


EST12 


a11. 12 


EST13 


b1,2 


CABI 


b3, 4 


EST17 


b5.6 


GA4 


b7,8 


EST19 


b9, 10 


GBF-1 


b11. 12 


EST23 


cl.2 


EST29 


C3.4 


GBF-2 


c5,6 


EST34 


c7, 8 


EST35 


c9. 10 


EST41 


c11, 12 


rGR 


d1 t 2 


EST42 


d3, 4 


EST45 


d5. 6 


HAT1 


d7 t 8 


EST46 


d9, 10 


EST49 


dl1, 12 


HAT2 


e1.2 


HAT 4 


e3. 4 


EST50 


e5. 6 


HATS 


e7, 8 


EST51 


e9, 10 


HAT22 


ell, 12 


EST52 


f1.2 


EST59 


f3.4 


KNAT1 


15.6 


EST60 


f7. 8 


EST69 


f9. 10 


PPH1 


f11, 12 


EST 70 


01.2 


EST75 


g3, 4 


EST 78 


g5,6 


ROC1 


g7,e 


EST82 


g9, 10 


EST83 


gn. 12 


EST84 


hl,2 


EST91 


h3, 4 


EST96 


h5,6 


SARI 


h7,8 


EST100 


n9, 10 


EST103 


h11, 12 


TRP4 



Function 



Human AChR 
Actin 

NADH dehydrogenase 

Actin 1 

Unknown 

Actin 

Chlorophyll a/b binding 
Phosphogtycerate kinase 
Gtbbereiiic acid biosynthesis 
Unknown 

G-box binding factor 1 
Elongation factor 
Aldolase 

G-box binding factor 2 
Chforopiast protease 
Unknown 
Cataiase 

Rat glucocorticoid receptor 

Unknown 

ATPase 

Homeobox-teucine zipper 1 
Light harvesting complex 
Unknown 

Homeobox-teucine zipper 2 
Homeobox-teucine zipper 4 
Phosphoribulokinase 
Homeobox-teucine zipper 5 
Unknown 

Homeobox-teucine zipper 22 
Oxygen evorving 
Unknown 

KnoffecMike horneobox 1 
RuBtsCO small submit 
Translation elongation factor 
Protein phosphatase 1 
Unknown 

Chloroptast protease 

Unknown 

Cyclophilin 

GTP binding 

Unknown 

Unknown 

Unknown 

Unknown 

Syrtaptobrevin 

Light harvesting complex 

Light harvesting complex 

Yeast tryptophan biosynthesis 



Accession 
number 



H36236 

227010 

M20016 

U36594t 

T45783 

M85150 

T44490 

L37126 

U36595t 

X63894 

X52256 

T04477 

X63895 

R87034 

T14152 

T22720 

M 14053 

U36596t 

J04185 

U09332 

T04063 

T76267 

U09335 

M90394 

T04344 

M90416 

233675 

U09336 

T21749 

234607 

U14174 

X14564 

T42799 

U34803 

T44621 

T43698 

R65481 

L14844 

X59152 

233795 

T45278 

T13832 

R64816 

M90418 

218205 

X03909 

X04273 



observed for HAT4, large changes in ex- 
pression were not observed for any of the 
other 44 genes we examined. This was 
somewhat surprising, particularly because 
comparative analysis of leaf and root tissue 
identified 27 differentially expressed genes. 
Analysis of an expanded set of genes may be 
required to identify genes whose expression 
changes upon HAT4 overexpression; alter- 
natively, a comparison of mRNA popula- 
tions from specific tissues of wild-type and 
HAT4-transgenic plants may allow identi- 
fication of downstream genes. 

At the current density of robotic printing, 
it is feasible to scale up the fabrication pro- 
cess to produce arrays containing 20,000 
cDNA targets. At this density, a single array 
would be sufficient to provide gene-specific 
targets encompassing nearly the entire rep- 
ertoire of expressed genes in the Arabidopw 
genome (2). The availability of 20,274 ESTs 
from Arabidopsis (1,9) would provide a rich 
source of templates for such studies. 

The estimated 100,000 genes in the hu- 
man genome (10) exceeds the number of 
Arabidopsis genes by a factor of 5 (2). This 
modest increase in complexity suggests that 
similar cDNA microarrays, prepared from 
the rapidly growing repertoire of human 
ESTs (I), could be used to determine the 
expression patterns of tens of thousands of 
human genes in diverse cell types. Coupling 
an amplification strategy to the reverse 
transcription reaction (J I) could make it 
feasible to monitor expression even in 
minute tissue samples. A wide variety of 
acute and chronic physiological and patho- 
logical conditions might lead to character- 
istic changes in the patterns of gene expres- 
sion in peripheral blood cells or other easily 
sampled tissues. In concert with cDNA mi- 
croarrays for monitoring complex expres- 
sion patterns, these tissues might therefore 
serve as sensitive in vivo sensors for clinical 
diagnosis. Microarrays of cDN As could thus 
provide a useful link between human gene 
sequences and clinical medicine. 



Table 2. Gene expression rnonitoring by rnicroar- 
ray and RNA blot analyses; tg, H4r4-transgenic. 
See Table 1 for additional gene information. Ex- 
pression levels (w/w) were calibrated with the use 
of known amounts of human AChR mRNA. Values 
for the rnicroarray were determined from rnicroar- 
ray scans Fig. 1); values for the RNA Wot were 
determined from RNA blots (Rg. 2). 



•Proprietary sequence of Stratagene (La Jotta, California). tNo match in the database; novel EST. 



Gene 


Expression level (w/w) 


Microarray 


RNA blot 


CABI 


1:48 


1:83 


CABt (tg) 


1:120 


1:150 


HAT 4 


1:8300 


1:6300 


HAT4 (tg) 


1:150 


1:210 


ROC1 


1:1200 


1:1800 


ROC1 (tg) 


1:260 


1:1300 
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Gene Therapy in Peripheral Blood 
Lymphocytes and Bone Marrow for 
ADA" Immunodeficient Patients 

Claudio Bordignon,* Luigi D. Notarangelo, Nadia Nobili, 
Giuliana Ferrari, Giulia Casorati, Paola Panina, Evelina Mazzolari, 
Daniela Maggioni, Claudia Rossi, Paolo Servida, 
Alberto G. Ugazio, Fulvio Mavilio 

Adenosine deaminase (ADA) deficiency results in severe combined immunodeficiency, 
the first genetic disorder treated by gene therapy. Two different retroviral vectors were 
used to transfer ex vivo the human ADA minigene into bone marrow cells and peripheral 
blood lymphocytes from two patients undergoing exogenous enzyme replacement ther- 
apy. After 2 years of treatment long-term survival of T and B lymphocytes, marrow cells, 
and granulocytes expressing the transferred ADA gene was demonstrated and resulted 
in normalization of the immune repertoire and restoration of cellular and humoral immunity. 
After discontinuation of treatment, T lymphocytes, derived from transduced peripheral 
blood lymphocytes, were progressively replaced by marrow-derived T cells in both pa- 
tients. These results indicate successful gene transfer into long-lasting progenitor cells, 
producing a functional multilineage progeny. 



Severe combined immunodeficiency asso- 
ciated with inherited deficiency of ADA 
( J ) is usually fatal unless affected children 
are kept in protective isolation or the im- 
mune system is reconstituted by bone mar- 
row transplantation from a human leuko- 
cyte antigen (HLA)-identical sibling donor 
(2). This is the therapy of choice, although 
it is available only for a minority of patients. 
In recent years, other forms of therapy have 
been developed, including transplants from 
haploidentical donors (3,4), exogenous en- 
zyme replacement (5), and somatic-cell 
gene therapy (6-9). 

We previously reported a preclinical mod- 
el in which ADA gene transfer and expression 
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successfully restored immune functions in hu- 
man ADA-def icient (ADA" ) peripheral 
blood lymphocytes (PBLs) in immunodefi- 
cient mice in vivo (10, 11). On the basis of 
these preclinical results, the clinical applica- 
tion of gene therapy for the treatment of 
ADA" SC1D (severe combined immunodefi- 
ciency disease) patients who previously failed 
exogenous enzyme replacement therapy was 
approved by our Institutional Ethical Com- 
mittees and by the Italian National Commit- 
tee for Bioethics (12). In addition to evaluat- 
ing the safety and efficacy of the gene therapy 
procedure, the aim of the study was to define 
the relative role of PBLs and hematopoietic 
stem cells in the long-term reconstitution of 
immune functions after retroviral vector-me- 
diated ADA gene transfer. For this purpose, 
two structurally identical vectors expressing 
the human ADA complementary DNA 
(cDNA), distinguishable by the presence of 
alternative restriction sites in a nonfunctional 
region of the viral long-terminal repeat 
(LTR), were used to transduce PBLs and bone 
marrow (BM) cells independently. This pro- 
cedure allowed identification of the origin of 
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